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Abstract. We study the formation of the spectra in ad- 
vection dominated accretion flows (ADAFs) around Kerr 
black holes. We use a Monte Carlo approach and fully 
general relativistic treatment to follow the paths of indi- 
vidual photons and model their scattering with mildly rel- 
ativistic, thermal electrons of the two temperature plasma 
present in the flow. We are mostly interested in the depen- 
dence of the spectra on the black hole angular momentum, 
and we find that the influence of the black hole rotation 
rate on the flow structure has an impact on the result- 
ing spectra. The flows around the fast rotating holes pro- 
duce relatively harder spectra. This property of the models 
should be taken into account when modeling the individ- 
ual sources and the population of inefficiently accreting 
black holes in the Universe. 
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1. Introduction 

The Advection Dominated Accretion Flows (ADAFs) 
have been reviewed by a number of authors, most recently 
by Narayan, Mahad evan, & Quat aert (1998 ), Kato, Fukue, 
& Mineshige ( |1998D , and Lasota ( |1998[ ). The ADAFs rep- 
resent a class of optically thin solutions of the accretion 
flow, which radiate so inefficiently, that almost all the heat 
dissipated inside the fluid is subsequently advected to- 
ward the black hole horizon. Since the cooling of matter 
is negligible, the equations of fluid dynamics are indepen- 
dent of the equations describing the emission, absorption 
and scattering of radiation. This allows one to address 
these two topics separately. In this paper we are concerned 
mostly with the radiation processes inside the flow. 

The fully relativistic dynamics of ADAFs has been 
described and some solutions have been obtained in La- 
sota ( |1994D , Abramowicz et al. ( |1996|, h ereafter ACGL), 
Abramowicz, Lanza & Percival ( |I997 ), Peitz & Appl 
( |I997D , Jaroszynski & Kur piewsk i ( 1997] , hereafter Pa- 
per I ), Ga mmie & Popham (1998), and Popham & Gam- 
mie (|l998). The most extensive survey of the parameter 



space is probably presented by Popham & Gammie ( |1998D , 
where the dependence of the flow on the black hole spin 
a, the viscosity parameter a, the gas adiabatic index and 
on the advection parameter / (where / = 1 means fully 
advective flow without any cooling) is investigated. This 
work shows rather strong dependence of the flow charac- 
teristics on parameters mentioned, especially on the black 
hole spin. The dependence on advection parameter is also 
substantial, but for a narrow range of /, which can rep- 
resent flows with negligible cooling (0.9 < / < 1) one 
can use the f — 1 solutions. Similarly, the solutions de- 
pend strongly on the viscosity parameter a, but for the 
limited range of this parameter (a > 0.01), which is more 
physically relevant (Balbus, Gammie & Hawley 1995| ), the 
differences between the solutions are not dramatic, in par- 
ticular the topology of the isobars is the same. 

In our calculations we use the solutions of the equa- 
tions describing the flow dynamics presented in Paper I. 
All our models use / = 1 and a = 0.1. Both values are rep- 
resentative of the physically relevant ADAFs. The black 
hole spin, we consider, is limited to the three values (a = 0, 
0.5, and 0.9). 

The main purpose of this paper is a self-consistent 
treatment of photon Comptonization in a two temperature 
plasma of an ADAF solution. To do so we need a 3D distri- 
bution of matter density, velocity and temperature, while 
the standard solutions give only the vertically averaged 
quantities as measured at the equator. At this point it is 
important to choose the relevant method of ve rtical aver- 
aging and w e follo w Abramowicz et al. ( 1997 ), Quataert 
& Narayan ( |1998|) and Paper I in choosing averaging on 
spheres (and not on cylinders) . In this approach the mat- 
ter distribution in space is limited by the centrifugal forces 
barrier and does not resemble the infinite isothermal at- 
mosphere, obtained in the alternative approach. 

We calculate the spectrum of photons leaving the flow. 
The photons originate in bremsstrahlung and thermal syn- 
chrotron processes, which can be described locally. The 
scattering (if any) can take place at any point of the pho- 
ton trajectory inside matter distribution and is nonlocal. 
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We simulate this process using Monte Carlo approach. 
Our simulation of photon Comptonization is quite stan- 
dard, but the flow, where the processes take place is rather 
complicated. Since the optical depth for scattering is low, 
a photon can travel to a distant part of the fluid before un- 
dergoing any interaction. The relative motion of the fluid 
elements, where consecutive interactions of photons with 
matter take place, can be substantial. Also, the photons 
traveling in the densest parts of the flow, near the horizon, 
can be deflected by the gravitational field. This and other 
relativistic effects in photon motion are included in our 
study. 

The main aim of our investigation is the self-consistent 
treatment of Comptonization. Most non-analytical ther- 
mal Comptonization models use iterative methods of solv- 
ing the kinetic equations (e.g. Poutanen & Svensson f 996| ). 
This method has been adopted by Narayan, Barret & Mc- 



Clintock (1997) to calculate the spectrum of soft X-ray 
transient source V404 Cyg in which, they believe, ADAF 
exists. In our calculations we use Monte Carlo method 
which has been described in most detail by Pozdnyakov, 



Sobol & Sunyaev ( [l977[ ) and Gorecki & Wilczewski (|I98J. 
This method allows us to follow only one photon at a 
time and we are not able to take into account the creation 
of e^e~ pairs. However as Bjornsson et al. ( 1996D and 
Kusunose & Mineshige ( |1996[ ) showed, the role of e+e~ 
pairs in ADAFs is not significant. 

In the next Section we briefly characterize the model 
of ADAF. The method of Monte Carlo Comptonization 
is presented in Sec. 3. In Sec. 4 we present the results of 
calculations, showing the ADAFs spectra. The discussion 
and conclusions follow in the last Section. 



2. The 3D description of the fluid 



2.1. The accretion flow 

We investigate the stationary flow of matter and the prop- 
agation of light in the gravitational field of a rotating Kerr 
black hole using the Boyer-Lindquist coordinates t, 4>, r, 6 
and the metric components gabifiO) as given by Bardeen 
( 1973 ). We follow the (—,+,+,+) signature convention. 
We use geometrical units, so the speed of light c = I 
and the mass of the hole M = 1. The Kerr parameter a 
(0 < a < 1) gives the black hole angular momentum in 
geometrical units. We use the Einstein summation con- 
vention, where needed, and a semicolon for the covari- 
ant derivative. The normalization of four velocity with the 
chosen metric signature reads u°'Ua = — 1. 

The system of equations we use follows in general 
ACGL. Since we are neglecting the cooling processes at 
this stage, and treat the accretion flow in the disk ap- 
proximation, two components of velocity and the speed of 
sound cs given as functions of radius, fully describe the 
dynamics. We use the angular velocity Q. and the physi- 
cal radial velocity V as measured by locally nonrotating 



observers as main kinematic variables (compare ACGL). 
After the vertical averaging the velocity perpendicular to 
the equatorial plane is neglected {u^ = 0) and other com- 
ponents of the four velocity are given as: 



1 



f V 



(1) 

(2) 
(3) 



The above velocity components have to be known only 
at the equatorial plane and its close vicinity if one needs to 
obtain a system of equations describing an ADAF in the 
slim disk approximation (Abramowicz et al. |l98§| ). To de- 
scribe the matter distribution in space (also far from the 
equatorial plane) one has to introduce further specifica- 
tions of the velocity and angular momentum distribution, 
since only the averaged values enter the equations. We as- 
sume that the poloidal velocity component (u^) is absent 
and that the radial velocity component V depends on the 
radius only: 



V = V{r) EE 



(4) 



(The BL velocity component depends on Q through g^r 
- compare Eq.3) 

The choice of angular momentum distribution is less 
obvious. Close to the horizon, where the velocities are 
highest and possible kinematic effects most important, the 
specific angular momentum i = —u^/ut is approximately 
constant (Paper I, ACGL, Peitz & Appl, |I997| ). We as- 
sume I to be exactly constant on spheres: 



£ = l{r) 



(5) 



Since the metric components do depend on 9 the angular 
velocity is not constant on spheres: 



g^'ut + g^^u^ _ g 



4>t 



g^*^ut + g^'I'u^ - 
The other kinematic quantities are given as: 
I 

ut = 



(6) 

(7) 
(8) 



We assume that the accreting plasma contains small 
scale, isotropically tangled magnetic field resulting from 
magnetohydrodynamical instability (Balbus et al. 
Hence we write the total pressure as 



1995) 



P=Pg+ Pm Pg = PP Pn 



Sir 



(1 - /3)p 



(9) 



where pg is the gas pressure, - the magnetic pressure, 
B - the magnetic field, and /3 is a constant parameter. The 
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pressure, the rest mass density po and the sound velocity 
cs are related: 



(10) 



In our calculations we use a two temperature plasma with 
a small amount of magnetic field to represent the matter 
properties. Thus the ion pressure dominates, the gas is 
non-relativistic and the specific enthalpy /i is given as: 



<^+P 1,52 

Po ^ 



(11) 



where e is the total (rest mass plus thermal) energy den- 
sity. 

In the ADAF set of equations only the vertically aver- 
aged sound speed is used. It is in spirit of our approxima- 
tions to postulate: 



Cs = cs(r) 



(12) 



which means that the gas is isothermal on spheres. The 
electron temperature is much lower than the ion tem- 
perature, so its influence on the equation of state and the 
structure of the flow can be neglected. 

Our detailed (but approximate) description of the fluid 
kinematics makes it possible to find the density depen- 
dence on the angular coordinate 9. The set of equa- 
tions describing ADAF contains viscosity terms, which 
are included in the energy equation, but neglected in the 
mechanical equilibrium equations (ACGL, Peitz & Appl 



1997] , Paper I). Thus it is sufficient to use the ideal fluid 



energy momentum tensor: 

T,^ ^ (e + p)u''ub + pS^ 

The 6 conservation equation, Tg.^ = reads: 



u Ue-a 



-P.fi 



e + p 

After some algebra we obtain: 

1 I V'^ 

{\nut)fi - ('"(''^ + 0^ cos^ S))fi 



(13) 



(14) 



1 - V 



2i-y2 



— ^(Inpo).e 



Since V does not depend on 9, the LHS of the equation is a 
full gradient of a quantity which can be called a potential 
tjj. This implies the solution for the density: 



Po(f) = Po,cqexp 



(16) 



where the po.oq, V'cq denote the values measured at the 
equatorial plane. 

The constant sound speed on the spheres may sug- 
gest that the exponential atmosphere never ends. For the 
rotating configuration there is, however, an infinite poten- 
tial barrier close to the rotation axis, where ut oo and 



Pq ~f 0. Thus the vicinity of the rotation axis is empty and 
the density falls steeply down near this region. In this re- 
spect the ADAF solutions are similar to the so called thick 



accretion disks (Abramowicz, Jaroszynski, & Sikora 1978). 
Both have empty funnels around their rotation axes. 

The mass flow through the r = const surface can be 
calculated as: 



M = - I d9d(f>^/^po{0)u'' 
la 



(17) 



= '^''\Iy~^ 'iS^Ji9tci>~ 9u9,i>4,) gee Paid) (18) 

Combining the last equation and the formula for the 9 
dependence of the density we obtain the equatorial value 
of the density. 



2.2. Two temperature plasma 

While modeling the dynamics of ADAF we neglect the 
heat transfer and all the radiation processes, assuming 
that only a small part of the total energy generated by 
viscous processes can be affected by them. Now we are 
going to model the radiation processes. 

We assume that the whole energy dissipated is trans- 
ferred to the ions. The ions heat the electrons via Coulomb 
collisions and the electrons lose their energy by the syn- 
chrotron, bremsstrahlung and inverse Compton cooling 
processes. In the ADAFs thermalization time-scale greatly 
exceeds the dynamical time-scale and the plasma remains 
two temperature. We find the electron and ion tempera- 
tures, Tc(r, 9) and Ti(r, 9), self-consistently using the equa- 
tion of state 



f3p^ 



pkTi pkTc 



(19) 



where pi = 1.29 and p^ = 1.18 arc effective molecular 
weights of the ions and electrons, and the condition of 
thermal equilibrium applied locally 



(20) 



(15) where is the rate of Coulomb heating of electrons 



by ions (e.g. Mahadevan 1997 ), and g^^. are the syn- 
chrotron and bremsstrahlung cooling rates, and q^(j and 
%T c Compton cooling rate of synchrotron and 

bremsstrahlung photons, respectively. 

The calculation of the synchrotron cooling rate is 
somewhat complicated because the optical depth to ab- 
sorption (Self Synchrotron Absorption) for majority of the 
synchrotron photons is high. In our approach we assume 
that a photon can either escape from the medium carry- 
ing out its energy and taking part in cooling process, or 
be absorbed very close to the point of its original emission 
and not contributing to the cooling. In reality all photons 
carry energy and some are absorbed far from the emission 



4 



A. Kurpiewski & M. Jaroszynski: Comptonization of photons in advection dominated accretion flows 



point transferring energy to distant portions of the fluid. 
Our assumption neglects the heat transport between dif- 
ferent volume elements, but makes calculations doable. 

To find the probability of a photon escape from given 
location we follow N ~ 10^ rays with directions randomly 
chosen at the frame comoving with the fluid. The optical 
depth along a ray measured at the frequency v is 



4:7: Bt^ (i^) 



dl 



(21) 



where eg is the synchrotron emissivity, Bt^ is the Planck 
function corresponding to the electron temperature Tc, 
and the integration is over the proper distance, Zoo cor- 
responding to the point on the fluid boundary. Since the 
integrand in the above formula is the function of the elec- 
tron temperature we must assume the approximate dis- 
tribution of the electron temperature in the ADAF. As a 
first approximation we use our results from Paper I based 
on the approach of Narayan & Yi 1994 . Averaging one 
gets the probability of escape: 



1 



e " = 



N 

N ^ 

i=l 



(22) 



Since we assume that cooling is provided only by the es- 
caping photons, we have: 



eg(i^) e '^'^ dv 



(23) 



for the synchrotron cooling. We adopt the expressions for 
synchrotron emissivity from Pacholczyk ( |1970| ) and Ma- 
hadevan, Narayan & Yi ( 1996| ). 

For the bremsstrahlung cooling the solution is straight- 
forward. The absorption of low frequency bremsstrahlung 
photons has no practical meaning, so we adopt their fre- 
quency integrated emission as cooling rate 



e^,{v) dv 



(24) 



We take the expression for bremsstrahlung cooling rate 



from Stepney & Guilbert (1983) 



Finally we take the cooling by Comptonization of both 
synchrotron and bre msstra hlung photons using the for- 
mulae of Esin et al. ( 1996 ). The mean optical depth to 
Compton scattering is calculated in the same manner as 
above. Solving the equations |l^ and we get the electron 
temperature and the ion temperature. 

2.3. Spectra neglecting Comptonization 

As a by-product of the calculations of the previous sub- 
sections one can obtain the spectrum of the model, which 
would be valid if the Comptonization were unimportant. 
For sufficiently low frequencies a photon is rather absorbed 
than scattered, so both approaches, neglecting and includ- 
ing Comptonization, should give similar results in this 



regime. This gives a chance of a self check of the simu- 
lations. 

Calculation of many rays sent from a given fluid ele- 
ment can also be used to find the contribution of this el- 
ement to the total luminosity of the configuration as seen 
by a distant observer. If the frequency Vem and the direc- 
tion of a photon in the fluid frame is known, its frequency 
in the Boyer-Lindquist coordinate frame I'obs can be cal- 
culated and this is the frequency that would be measured 
by a distant observer, unless the photon goes under the 
horizon. For photons going to infinity the redshift factor 
can be defined: 



l + z- 



(25) 



Distant observers measure photon energies divided by the 
factor (l + z). Also the time interval between the detec- 
tion of two signals is the interval between their sending 
multiplied by (1 -I- z). Thus the contribution to the to- 
tal luminosity from the fluid element of the volume AV^, 
measured by distant observers in their frequency interval 



bs IS given as: 



i(i/obs)Az^obs ^ JT 2^ ^ AV" (26) 



N 

i=l 

where for the i-th ray: 



I'cm = i^obs(l + Z,) Al/cm = Aj/obs(l + Z,) 



(27) 



The summation is limited to rays which reach infinity. 
Integration over the volume of the fluid gives the total 
luminosity of the model. 

The emission from the configuration is not isotropic, 
so the observers at different position angles measure a dif- 
ferent flux of radiation. The luminosity calculated above 
is in fact an average of luminosities assigned to the disk by 
observers uniformly distributed on a sphere around the ob- 
ject. One can also find the average luminosity that would 
be measured by observers from a limited solid angle Af2. 
To do so it is sufficient to neglect in the summation all the 
rays which do not enter the region of interest and multiply 
the result by the correction factor Aw/ Aft. 

2.4- Spatial distribution of the photon emissivity 

In the Monte Carlo simulations we follow individual pho- 
tons as they travel through the fiuid undergoing consecu- 
tive scatterings. We have to know what is the distribution 
of the points of emission of the photons. We divide the flow 
into several spherical layers. The radius in the middle of 
the layer numbered j is rj . Each layer is then subdivided 
into annuli of limited range in the polar angle 9 between 
fluid boundaries Ominirj) and Tr—Ominifj). The angular co- 
ordinate at the middle of the annulus numbered jk is 9jk- 
One can assume that all fluid parameters are almost uni- 
form inside each annulus, and take their values at (rj, Oj^.) 
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as representative. Using similar arguments as in the pre- 
vious subsection we calculate the rate of photon emission 
from the region numbered jk of the volume AVjk- 



jk 



N 



N 

E 

1=1 



.(1 



(28) 



We adopt the expressions for synchrotron cooling from 
Pacholczyk ( |1970D and Mahadevan, Narayan fc Yi (|1996D , 
and for bremsstrahlung cooling from Svensson ( |1982| ) [and 
references therein]. The cooling rates are functions of the 
electron temperature, the number density of ions and elec- 
trons and the magnetic field density (synchrotron radia- 
tion) and hence they are functions of r and 9. We take 
e = e{rj,9jk)- The expression under the integral is regu- 
lar at low frequencies despite the presence of I'em in the 



denominator because Tj 



(i) 



oo when j^e 



0. The red- 



shift factor in the denominator takes care of the difference 
between clock rates in the fluid frame and at infinity. 

3. The Comptonization 

We follow the method of Comptonization described by 



Gorecki & Wilczewski (1984) 



3.1. Basic concepts 

The differential cross section for Compton scattering is 
given by the following formula (Akhiezer & Berestetski 



1965) 



— = 4^(l--0/c)-f— V 



(29) 



where hv, r2, hv' and Q,' are respectively the energy and 
the direction of the photon before and after the scattering, 
V is the velocity of the electron, 7 is the Lorentz factor and 
rg is the classical electron radius. The symbol X denotes 
the invariant part of the cross section : 



X 



X 



X' 
X 



where 

X hv , „ , X 
- = — 7(1 - t;J^/c) 
2 mc 



X' 

2 



hv' 



mc 



.7(1 - vn'/c) 



(30) 



(31) 



are the energies of the incoming photon and of the scat- 
tered photon, respectively, expressed in units of mc^ in 
the reference frame of the electron. The energies hv and 
hv' are related by the Compton formula 

hv{l - vn/c) 



hv' 



-^(l-Of2') 



(32) 



1 - vQ.' jc- 

We use the total Compton cross s ection which is given by 
(Berestetski, Lifshitz & Pitaevski 1972 ) 

a{x) = 



2^r2i 



(33) 



The basic concept of this method is to follow the photon 
trajectory from the moment of emission until the photon 
leaves the flow. The probability that a photon leaves the 
flow without scattering is 



Pi — exp 



ne{a)dl 



(34) 



where the integral is taken along the photon trajectory 
from the point of the last (i-th) scattering to the boundary 
of the flow Too, ne = J ne{v)d^v is the electron density. 



ne{v) is the electron velocity distribution, and 

(a) = — [ nJv){l - v^l/c)a(x)d^ 
ne J 



(35) 



is a mean cross section averaged over the electron velocity 
distribution. ne{v). 

The probability Pi enables us to flnd the statistical 
weights of the number of photons leaving the flow without 
scattering (and thus contributing to the emerging spec- 
trum) and the photons which remain in the flow and un- 
dergo the next {[i + l]-th) scattering. These weights are 
given by WiPi and Wi+i = Wi{l — Pi), respectively, where 
i = 0, 1, 2, 3, ... is the index denoting the succeeding 
scatterings. We assume wq = 1. We follow the trajectory 
of the photon until w becomes less than a certain mini- 
mal value Wniin- Since ADAFs are optically thin (in our 
model the Thomson optical depth is about 0.1 in equato- 
rial directions) the mean number of scatterings is 4 - 5 for 
Wmin = 10""^. 

3.2. Generating the random variables 

The random variables are generated from the probability 
distributions using Monte Carlo methods: the inversion of 
the cumulative distribution function or von Neumann's 
rejection technique. Multi-dimensional distributions are 
modeled using the conditional probability distributions. 

(i) At flrst we generate the position vector at which the 
photon is initially emitted. According to our approxima- 
tion this position is uniformly distributed in space within 
each of the annuli jk. The probability that a photon is 
emitted from a region of the given number jfc is : 



fjk 



jk 



T^k'^r^ 



(36) 



where jifjk are given in Eq. In Eq. |2^ we take into 
account the optical depth due to the absorption, so the 
number of the low frequency photons we use in the simu- 
lation is the expected number of the photons which have 
a chance of escape. The absorption does not have to be 
considered on the further photon trajectory. (After a scat- 
tering with relativistic electrons a photon gains so much 
energy that the possibility of its absorption can be ne- 
glected. The probability of absorption on the original tra- 
jectory is included in Eq. p8| .) 
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(ii) As the position of input photon is determined we 
can generate the initial energy of the photon hvQ using a 
probability distribution specified by the photon spectrum 
of synchrotron or bremsstrahlung emission 



(37) 



is the photon spectrum 



where n„{rj,9jk)[Hz^^cm^^ 
approximated by formulae of Pacholczyk (1970) and Ma- 
hadevan, Narayan & Yi ( 1996| ) for synchrotron emission 
or Svensson (1982) for bremsstrahlung emission. The pho- 
ton spectrum is determined from the energy spectrum by 
dividing the last one by hv. 

(iii) We assume that the emission of input photons is 
isotropic. Hence the direction of the photon in a comoving 
Cartesian coordinate frame fl— (sin 8 cos $, sin0sin$, 
cos Q) is generated from uniform distributions in the 
ranges cos 6 e [-1; 1] and $ e [0; 27r]. 

In this way we determine the set of parameters {tq, 
hvo, f2o, wq = 1} describing the initial point of the trajec- 
tory of the photons beam. We calculate the next points of 
the trajectory, i.e. {r^+i, hv^+i, $7^+1, w^+i} {i = 0, 1, 2, 
3, ...) until w = Wmin- Thc way of computing the weights 
Wi is described in Sec. 3.1. Below we present the way of 
computing the position, energy and direction of a photon 
after following scatterings. 

(iv) The position ri+i is found on thc photon trajec- 
tory at the proper distance I from the starting point r^, 
from the probability distribution: 



where 



,-r(n dr(;) 
dl 



dl 



t{1) = / n4c7)dl 
Jq 



(38) 



(39) 



(v) The two remaining parameters hvi+i and ili+i of 
the {i + l) — th point of the photon trajectory are obtained 
by simulating the scattering of the photon of energy hvi 
and direction $7^ by an electron with velocity v. To de- 
scribe the probability distribution of this scattering we 
use the differential cross section (1) : 



n{v){\ — vVli/c) 



do- 



//n(t;)(l - vrL,/c)^d^n'd?v 



(40) 



We model the multi-dimensional probability distribution 
(13) as a product of the probability distribution of v and 
the conditional probability distribution of ^i+i 



where 

hi.v) = — (1 - c)—— 

Tie W 



(41) 



(42) 



and 



\v) = 



1 



da 



a{x) dQi+i 

We generate the velocity v from Eq 
from Eq 



(43) 



direction Cli 



43 



3 and then the 
Having Jli+i , the energy 



of the scattered photon hvi+i can be obtained from the 
Compton formula (|3^). The detailed description of the 
method of modeling the probability distributions fi(v) 
and f 2{^i+i\v) can be found in Gorecki & Wilczewski 
( |1984D . 



4. Results 

We have performed trial calculations of the ADAFs spec- 
tra employing thc method described. We use the models of 
ADAFs from Paper I. For the black hole mass, accretion 
rate and parameter P we use the parameters of Lasota et 
al. (|1996| ), Mbh = 3.6 x IO^Mq, m = 0.016, l3 = 0.95, 
which they apply in the modeling of NGC 4258. We do 
not, however, attach cold, thin disks at large radii to our 
ADAF solutions. In Fig. |] we show the input spectra of 
synchrotron and bremsstrahlung photons as well as the re- 
sulting Comptonized spectrum of the disk around a = 0.9 
black hole. All the spectra in this and other diagrams are 
shown as \g{v) versus \g{i'F^) plots. For the synchrotron 
input the results are based on calculations including 10^ 
input photons and following more than 5 x 10^ branches 
of photon trajectories. The Comptonization plays a less 
important role in the case of bremsstrahlung radiation, so 
we use ~ 10 times fewer photon trajectories to obtain the 
spectra in this case. Since the Comptonization preserves 
the photon number, we normalize the spectra using Eq. |2^ 
with either synchrotron or bremsstrahlung emissivity un- 
der the integral to obtain the relative numbers of seed 
photons of each kind. 

The Comptonized spectra of synchrotron radiation are 
smooth enough to allow a power law fit with power index 
r (i.e. Li, ~ t^^^)- The fit is not valid at the vicinity of 
the first peak, which is due to the seed photons. In Fig. || 
we show synchrotron spectra with fits. As can be seen in 
the plots the slopes of the spectra depend on the model. 
Since the black hole mass and the accretion rate are the 
same for all three cases, the differences must be attributed 
to the black hole angular momentum and its influence on 
the flow structure. The power law indices estimated from 
fits are F = 0.89, 0.85, and 0.81 for a ^ 0., 0.5, and 0.9 
respectively. 

The three spectra resulting from combined effects of 
synchrotron and bremsstrahlung emission with Compto- 
nization are shown in Fig. ^. 

We have also checked the dependence of the observed 
total luminosity of our models on the observer's position. 
For the bremsstrahlung photons the dependence is absent. 
The synchrotron radiation observed from the equatorial 
plane is stronger by 10 to 20% as compared to the mea- 
surement from the axis of rotation. 
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Fig. 1. The spectra of the input synchrotron and brems- 
strahlung photons and the resuhing spectrum after Comp- 
tonization for a = 0.9. The resuhing spectrum is shown as 
a sohd Hue. Synchrotron input (left) and bremsstrahlung 
(right) use dotted hues. 



Following the photons we are able to find the frac- 
tion which goes under the horizon. For the synchrotron 
photons the numbers are: 0.07, 0.05, and 0.04 for a = 0, 
0.5 and 0.9 respectively. The bremsstrahlung photons are 
emitted at relatively larger distances from the horizon and 
less than 1% of them are lost in all cases. The fraction of 
photons emitted by the fluid and going under the horizon 
is a decreasing function of the black hole angular momen- 
tum according to our simulations. We have checked this 
result of our simulations making an independent calcula- 
tion. We have compared isotropic sources of radiation at 
the same distance from the black hole, comoving with the 
matter of the three models we use. Again the fraction of 
rays going under the horizon is the smallest for the case 
of the most rapidly rotating hole. The effect must be at- 
tributed to the differences in matter kinematics between 
the three models. 

5. Discussion and conclusions 

We use the ADAF models from Paper I, but we apply 
some changes to the treatment of the 3D structure of the 
flow. The models are still based on vertically averaged, 
stationary equations for the radial structure of the flow, 
which uses the vertical scale height and quantities mea- 
sured at the equator as flow variables. Such treatment can 
be inadequate near the flow boundary, especially close to 
the rotation axis. We assume the specific angular momen- 
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Fig. 2. The power law fits to the Comptonized syn- 
chrotron emission for models with black hole angular mo- 
mentum a — (dashed lines), a — 0.5 (dotted) and a = 0.9 
(solid). The corresponding straight lines represent the fits. 

tum of matter, its radial velocity, and the sound speed 
to be constant on spheres. Under such assumptions it is 
possible to introduce an effective potential on each sphere, 
which acts as infinite centrifugal potential barrier near the 
rotation axis. This infinite barrier causes the sharp drop 
to zero of matter density there despite the fact that the 
gas is approximately isothermal on spheres; it effectively 
removes matter with artificially high angular velocity. 

Our models are based on the assumption that all the 
heat dissipated in the flow goes into the ions. We ne- 
glect the direct viscous heating of electrons and the fact 
that their entropy changes as they fall toward the black 
hole, which represents the advection of heat. Both mecha- 
nisms influence the energy balance equation and, as shown 
by Nakamura et al. (1997), Narayan et al. (1998) and 
Quataert & Narayan (1998), have an impact on the elec- 
tron temperature and the resulting spectrum of the model. 
While advection of heat by electrons is a well defined pro- 
cess, the viscous heating must be introduced using another 
free parameter. The combined effects of viscosity and ad- 
vection on electrons would infiuence all our models in a 
similar way, not greatly changing the differences between 
them. 

In our calculations we use three models of the matter 
flow onto the black hole, with the same accretion rate and 
the same black hole mass, but with different black hole 
angular momentum. The spin of the black hole has strong 
influence on the density and temperature of the matter 
near the horizon, which are both increasing functions of 
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motion of the starting point of a photon and the plaee of 
its interaction with the electrons. 

We are not attempting to model NGC 4258, but using 
the parameters from the model of Lasota et al. (1996) 



we get the right luminosity in the X-rays for the a — 
0.9 model. The slope of the calculated spectrum at this 
frequency (10"'^* Hz) is within the observational bounds. 
The spectrum of NGC 4258 has been also modeled by 
Gammie, Narayan, & Blandford (199S). Although they 



use slightly different values of the ADAF paramet ers, th eir 
results are very similar to those of Lasota et al. ( 1996|) . 
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